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Abstract 

We measure thermodynamic magnetization of a low-disordered, strongly correlated 
two-dimensional electron system in silicon. Pauli spin susceptibility is observed to 
grow critically at low electron densities — behavior that is characteristic of the ex- 
istence of a phase transition. A new, parameter-free method is used to directly de- 
termine the spectrum characteristics (Lande (7-factor and the cyclotron mass) when 
the Fermi level lies outside the spectral gaps and the inter-level interactions between 
quasiparticles are avoided. It turns out that, unlike in the Stoner scenario, the crit- 
ical growth of the spin susceptibility originates from the dramatic enhancement of 
the effective mass, while the enhancement of the 5- factor is weak and practically 
independent of the electron density. 
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1 Introduction 



Presently, theoretical description of interacting electron systems is restricted 
to two limiting cases: (i) weak electron-electron interactions (small ratio of the 
Coulomb and Fermi energies = Ec/Ep <S 1, high electron densities) and 
(ii) very strong electron-electron interactions (r^, ^ 1, very low electron den- 
sities). In the first case, conventional Fermi-liquid behavior [1] is established, 
while in the second case, formation of the Wigner crystal is expected [2] (for 
recent developments, see Ref. [3]). Numerous experiments performed in both 
three- (3D) and two-dimensional (2D) electron systems at intermediate inter- 
action strengths (1 ^ ^ 5) have not demonstrated any significant change in 
properties compared to the weakly-interacting regime (see, e.g., Refs. [4,5]). It 
was not until recently that qualitative deviations from the weakly-interacting 
Fermi liquid behavior (in particular, the drastic increase of the effective elec- 
tron mass with decreasing electron density) have been found in strongly cor- 
related 2D electron systems (r^, > 10) [6]. However, these findings have been 
based solely on the studies of a kinetic parameter (conductivity), which, in 
general, is not a characteristic of a state of matter. 

The 2D electron system in silicon turns out to be a very convenient object for 
studies of the strongly correlated regime due to the large interaction strengths 
{ts > 10 can be easily reached) and high homogeneity of the samples esti- 
mated (from the width of the magnetocapacitance minima in perpendicular 
magnetic fields) at about 4 x 10^ cm~^ [7] . Here we report measurements of the 
thermodynamic magnetization and density of states in such a low-disordered, 
strongly correlated 2D electron system in silicon. We concentrate on the metal- 
lic regime where conductivity a ^ /h. We have found that in this system 
subjected to parallel magnetic fields, the spin susceptibility of band electrons 
(Pauli spin susceptibility) becomes enhanced by almost an order of magni- 
tude at low electron densities, growing critically near a certain critical density 
rZ;^ ~ 8 X 10^*^ cm~^: behavior that is characteristic in the close vicinity of a 
phase transition. The density is coincident within the experimental uncer- 
tainty with the critical density Uc for the zero-field metal-insulator transition 
(MIT) in our samples. We have also found by measurements in perpendicular 
and tilted magnetic fields that the (^-factor is weakly enhanced and practically 
independent of the electron density down to the lowest densities reached, while 
the cyclotron mass becomes strongly enhanced at low n^. Thus, unlike in the 
Stoner scenario, it is the effective mass, rather than the g-factor, that is re- 
sponsible for the dramatically enhanced spin susceptibility at low electron 
densities. 



2 Experimental setup and samples 

Magnetization is one of the least studied properties of 2D electron systems: 
signals associated with it are weak making measuring them a challenging ex- 
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periment. Few experimental observations of the de Haas- van Alphen effect in 
2D electron systems were made using SQUID magnetometers [8], pick up coils 
lithographed above the gate [9], or torque magnetometers [10]. Here we use 
a novel method [11] to measure the magnetization that entails modulating 
the magnetic field with an auxiliary coil and measuring the AC current in- 
duced between the gate and the 2D electron system. The experimental set-up 
is shown in Fig. 1 (a). Magnetic field B was modulated with a small ac field 
6B in the range 0.005 - 0.03 T at frequencies between / = 0.05 and 0.45 Hz. 
The in-phase and out-of-phase components of the current between the gate 
and the 2D electron system were measured with high precision (~ 10^^^ A) 
using a current-voltage converter and a lock-in amplifier. The real (in-phase) 
component of the current depends on the dissipative conductivity a of the 
sample. We are, however, interested in the imaginary (out-of-phase) current 
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Fig. 1. (a) Experimental set-up for measurements of thermodynamic magnetization, 
(b) Top view of a MOSFET with slit in the gate metalization indicated. 
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component, which, under condition that 27r/C <^ cr, is equal to 



. 27r/C d^M 

Im z = — 

e aiJ 



where C is the capacitance of the sample measured in the same experiment and 
II is the chemical potential. Use of ultra-low frequencies ensures that 2-71 fC <^ 
a and the out-of-phase current component is not contaminated by lateral 
transport effects. By applying the Maxwell relation dM/dus = —dfi/dB, one 
can then obtain the magnetization M from the measured i. For measurements 
of the thermodynamic density of states, a similar circuit was used with a dis- 
tinction that the gate voltage was modulated and thus the imaginary current 
component was proportional to the capacitance. Thermodynamic density of 
states dus/dii is related to magnetocapacitance via 



1 _ 1 1 

C ~ C~o^ Ae'^idnJdnY 



where Co is the geometric capacitance and A is the sample area. 

Measurements were made in an Oxford dilution refrigerator on low-disordered 
(lOO)-silicon metal-oxide-semiconductor field-effect transistors (MOSFETs). 
These samples are remarkable by the absence of a band tail of localized elec- 
trons down to electron densities Ug ^ 1 x 10^^ cm"^ [6], which allows one 
to study properties of a clean 2D electron system without admixture of local 
moments [12,13,14]. The level of disorder is determined by the quality of the 
Si-Si02 interface. The silicon oxide layer in our samples was grown repeat- 
edly during several stages of the fabrication to decrease number of charged 
impurities and obtain high quality interface; as a result, in samples with ox- 
ide thickness of 149 nm, we have achieved peak mobilities of 3 m^/Vs at 
T — 100 mK. The second advantage of these samples is a very low contact 
resistance (in "conventional" silicon samples, high contact resistance becomes 
the main experimental obstacle in the low density/low temperature limit). To 
minimize contact resistance, thin slits in the gate metalization have been in- 
troduced (see Fig. 1 (b)), which allows for maintaining high electron density 
near the contacts regardless of its value in the main part of the sample. 
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Fig. 2. Imaginary current component in the magnetization experiment as a function 
of the electron density in a magnetic field of 5 T and T = 0.4 K. Grey area depicts 
the insulating phase. Magnetization versus is displayed in the inset. Note that 
the maximum M is coincident within the experimental uncertainty with ^b^s- 

3 Experimental results and discussion 

3.1 Pauli spin susceptibility 

We start with measurements of Pauli spin susceptibility. A typical experimen- 
tal trace of i{ns) in a parallel magnetic field of 5 T is displayed in Fig. 2. 
The inset shows magnetization M{ns) in the metallic phase obtained by in- 
tegrating dM/dug = —dfj,/dB with the integration constant M{oo) = Bxq, 
where Xo is the Pauli spin susceptibility of non-interacting electrons. A nearly 
anti-symmetric jump of i{ns) about zero on the y-axis (marked by the black 
arrow) separates the high- and low-density regions in which the signal is pos- 
itive and negative (M(ns) is decreasing and increasing), respectively. Such 
a behavior is expected based on simple considerations. At low densities, all 
electrons are spin-polarized in a magnetic field, so for the simple case of non- 
interacting 2D electrons one expects dfi/dB = —jiB (at — > 0, deep in the 
insulating regime, the capacitance of the system vanishes and, therefore, the 
measured current approaches zero). At higher densities, when the electrons 
start to fill the upper spin subband, M(ns) starts to decrease, and dfj,/dB is 
determined by the renormalized Pauli spin susceptibility x is expected to 
decrease with due to reduction in the strength of electron-electron inter- 
actions. Finally, in the high-density limit, the spin susceptibility approaches 
its "non-interacting" value xo, and dfi/dB should approach zero. The onset of 
complete spin polarization — the electron density Hp at which the electrons 
start to fill the upper spin subband — is given by the condition dfi/dB = 
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Fig. 3. (a) The experimental d^/dB as a function of electron density in different 
magnetic fields and T = 0.4 K. The curves are vertically shifted for clarity. Grey 
area depicts the insulating phase. Note that the onset of full spin polarization in 
our experiment always takes place in the metallic regime, (b) Scaling of the d^/dB 
curves, normalized by magnetic field magnitude, at high electron densities. The 
dashed line represents the "master curve" . Spin susceptibility obtained by integrat- 
ing the master curve (dashed line) and raw data at B = 1.5, 3, and 6 T is displayed 
in the inset. 

{Mijis) reaches a maximum), as indicated by the black arrow in the figure. 
It is important that over the range of magnetic fields used in the experiment 
(1.5-7 tesla), the maximum M coincides within the experimental uncertainty 
with ^B^s thus confirming that all the electrons are indeed spin-polarized be- 
low Up. Note however that the absolute value of d/i/dB at Ug < Uc is reduced 
in the experiment. We attribute this to smearing of the minimum in i{ns) 
caused by possible influence of the residual disorder in the electron system, 
which is crucial in and just above the insulating phase, in contrast to the clean 
metallic regime we focus on here. Another reason for the reduction in dfi/dB 
is the electron-electron interactions (due to, e.g., the enhanced effective mass). 

In Fig. 3(a), we show a set of curves for the experimental dfi/dB versus elec- 
tron density in different magnetic fields. Experimental results in the range 
of magnetic fields studied do not depend, within the experimental noise, on 
temperature below 0.6 K (down to 0.15 K which was the lowest tempera- 
ture achieved in this experiment). The onset of full spin polarization shifts 
to higher electron densities with increasing magnetic field. Grey area depicts 
the insulating phase, which expands somewhat with B (for more on this, see 
Ref. [15]). Note that the range of magnetic fields used in our experiment is 
restricted from below by the condition that dfi/dB crosses zero in the metallic 
regime. In Fig. 3(b), we show how these curves, normalized by magnetic field, 
collapse in the partially-polarized regime onto a single "master curve". The 
existence of such scaling verifies proportionality of the magnetization to B, 
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Fig. 4. (a) Magnetocapacitance versus electron density for different magnetic fields, 
(b) Deviation of the C{ns) dependences for different magnetic fields from the B = 
reference curve. The traces are vertically shifted for clarity. The onset of full spin 
polarization is indicated by arrows. 

confirming that we deal with Pauli spin susceptibility of band electrons, and 
establishes a common zero level for the experimental traces. Integration of the 
master curve over yields the spin susceptibility x = ^/B, as shown in the 
inset to Fig. 3(b). Also shown is the spin susceptibility obtained by integration 
of raw curves at i? = 1.5, 3, and 6 tesla, which, within the experimental error, 
yield the same dependence. 



This method of measuring the spin susceptibility, being the most direct, suf- 
fers, however, from possible influence of the unknown diamagnetic contribution 
to the measured dii/dB, which arises from the finite width of the 2D electron 
layer [16]. To verify that this influence is negligible in our samples, we em- 
ploy another two independent methods to determine x- The second method is 
based on marking the electron density Hp at which dn/dB = and which cor- 
responds to the onset of complete spin polarization, as mentioned above. The 
so-determined polarization density np{B) can be easily converted into 
via X = fJ'B'np/B. Note that in contrast to the value of d/i/dB, the polarization 
density Up is practically not affected by possible influence of the diamagnetic 
shift. 



The third method for measuring Hp and insensitive to the diamagnetic shift, 
relies on analyzing the magnetocapacitance, C. Experimental traces C{ns) are 
shown in Fig. 4(a) for different magnetic fields. As the magnetic field is in- 
creased, a step-like feature emerges on the C{ns) curves and shifts to higher 
electron densities. This feature corresponds to the thermodynamic density of 
states abruptly changing when the electrons' spins become completely polar- 
ized. To see the step-like feature more clearly, in Fig. 4(b) we subtract the 
C{ns) curves for different magnetic fields from the reference B = curve. The 
fact that the jumps in C (as well as in dfi/dB) are washed out much stronger 
than it can be expected from possible inhomogeneities in the electron density 
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Fig. 5. Dependence of the Pauli spin susceptibility on electron density obtained 
by all three methods described in text: integral of the master curve (dashed line), 
dfi/dB = (circles), and magnetocapacitance (squares). The dotted line is a guide 
to the eye. Also shown by a solid line is the transport data of Ref. [7]. Inset: po- 
larization field as a function of the electron density determined from the magne- 
tization (circles) and magnetocapacitance (squares) data. The symbol size for the 
magnetization data reflects the experimental uncertainty, and the error bars for the 
magnetocapacitance data extend to the middle of the jump in C. The data for 
are described by a linear fit which extrapolates to a density close to the critical 
density for the = MIT. 

distribution (about 4 x 10^ cm~^ [7]) points to the importance of electron- 
electron interactions. Since the effects of interactions are different in the fuUy- 
and partially-polarized regimes, it is natural to mark the onset of full spin 
polarization at the beginning of the interaction-broadened jump, as indicated 
by arrows in the figure. In case the residual disorder does contribute to the 
jump broadening, we extend error bars to the middle of the jump, which yields 
an upper boundary for the onset of full spin polarization. 

In Fig. 5, we show the summary of the results for the Pauli spin suscepti- 
bility as a function of n<j, obtained using all three methods described above. 
The excellent agreement between the results obtained by all of the methods 
establishes that a possible influence of the diamagnetic shift is negligible [17] 
and, therefore, the validity of the data including those at the lowest electron 
densities is justified. There is also good agreement between these results and 
the data obtained by the transport experiments of Ref. [7]. This adds cred- 
ibility to the transport data and confirms that full spin polarization occurs 
at the field B^., however, we note again that evidence for the phase transition 
can only be obtained from thermodynamic measurements. The magnetization 
data extend to lower densities than the transport data, and larger values of x 
are reached, exceeding the "non-interacting" value xo by almost an order of 
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magnitude. The Pauli spin susceptibility behaves critically close to the critical 
density for the B = metal- insulator transition [18]: x ^s/ ("s — ""x)- This 
is in favor of the occurrence of a spontaneous spin polarization (either Wigner 
crystal [19] or ferromagnetic liquid) at low Ug, although in currently available 
samples, the formation of the band tail of localized electrons at Ug ^ ric con- 
ceals the origin of the low-density phase. In other words, so far, one can only 
reach an incipient transition to a new phase. 

The dependence Bc{ns), determined from the magnetization and magnetoca- 
pacitance data, is represented in the inset to Fig. 5. The two data sets coincide 
and are described well by a common linear fit which extrapolates to a density 
close to ric- We emphasize that in the low-field limit {B < 1.5 tesla), the 
jump in dfi/dB shifts to the insulating regime, which does not allow us to ap- 
proach closer vicinity of riy^: based on the data obtained in the regime of strong 
localization, one would not be able to make conclusions concerning properties 
of a clean metallic electron system which we are interested in here. Clearly, 
the fact that the linear Bc{ns) dependence persists down to the lowest electron 
densities achieved in the experiment confirms that we always deal with the 
clean metallic regime. 

In the end of this subsection, we would like to clarify the principal difference 
between our results and those of Ref. [11]. In the sample used by Prus et al, 
the critical density Uc for the B = MIT was considerably higher than in 
our samples caused by high level of disorder, and the band tail of localized 
electrons was present at all electron densities [11]. As a result, the crucial 
region of low electron densities, in which the critical behavior of the Pauli spin 
susceptibility occurs, falls within the insulating regime where the physics of 
local moments dominates [12,13,14]. Indeed, Prus et al. have found sub-linear 
M{B) dependence characteristic of local moments, and the extracted spin 
susceptibility in their sample has a Curie temperature dependence [14] . This 
is the case even at high electron densities, where metallic behavior might be 
expected instead. Such effects are absent in our samples: the spin susceptibility 
(in the partially-polarized system) is independent of the magnetic field and 
temperature, confirming that we deal with Pauli spin susceptibihty of band 
electrons. 



3.2 Magnetization in perpendicular magnetic fields: g -factor and effective 
mass 

Typical experimental traces of the gate current in a perpendicular magnetic 
field of 5 T are displayed in Fig. 6. Sharp dips in the out-of-phase component, 
seen at integer filling factors u = nshc/eB±, refiect gaps in the density of 
states: dips at odd filling factors correspond to the valley sphtting, the ones 
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Fig. 6. Out-of-phase (solid line) and in-phase (dotted line) current components 
as a function of the electron density in a perpendicular magnetic field of 5 T and 
T = 0.8 K. 6B = 0.022 T and / = 0.45 Hz. The value dfi/dB is indicated in units 
of the Bohr magneton //g. In the right-hand inset, we demonstrate proportionality 
of Im i to frequency: the solid and dashed lines (vertically shifted for clarity) corre- 
spond to 0.45 and 0.1 Hz, respectively; the y-component of the latter is multiplied 
by 4.5. The left-hand inset illustrates magnetization per electron. 

at z/ = 2 and 6 are due to the spin splitting, and the dip at u = 4 is due to 
the cyclotron splitting. However, there are no corresponding features in the in- 
phase current component, which ensures that we reach the low-frequency limit 
and the measured d/i/dB is not distorted by lateral transport effects. This is 
further confirmed by the fact that the out-of-phase current is proportional 
to the excitation frequency as displayed in the right-hand inset to Fig. 6. 
Magnetization per electron can be extracted by integrating the measured out- 
of-phase signal with respect to Ug, as shown in the left-hand inset to Fig. 6 for 
illustration. The magnetization exhibits the expected sawtooth oscillations, 
with sharp jumps at integer filling factors (note that the height of the jumps 
yields values that are smaller than the level splitting by the level width). 

If the disorder and interactions are disregarded, in quantizing magnetic fields 
(except at integer filling factors) the derivative dfi/dB = —dM/dng is equal 
to 



dB 



+ N 



2me 1 
± o^o 



(3) 



where /is is the Bohr magneton, iV is the Landau level number, me and 
rrih = 0.19 me are the free electron mass and band mass, respectively, and 
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Fig. 7. (a) Capacitance in B = 8 tesla and in S = as indicated. The (noise 
averaged) geometric capacitance is depicted by a dashed hne. (b) Im i oc d^/dB 
in a perpendicular magnetic field of 8 tesla. The maximum values possible in a 
non-interacting system (see text) are depicted by dashed lines. 



(7o = 2 is the gf-factor in bulk silicon. Disorder smears out the dependences 
which otherwise would consist of a series of delta-functions. Interactions mod- 
ify this picture in two ways: (i) by renormalizing the values of the cyclotron 
mass and ^^-factor and (ii) by providing a negative contribution of order 
— {e'^/elB){i^}^^'^ to the chemical potential [21,22] (here e is the dielectric con- 
stant, Ib is the magnetic length, and {v} is the deviation of the filling factor 
from the nearest integer). The latter effect, which is caused by the intra- level 
interactions between quasiparticles, leads to the so-called negative thermody- 
namic compressibility near integer filling factors predicted by Efros [22] and 
experimentally observed in Refs. [23,24]. 

In Fig. 7, we compare capacitance C with dfi/dB, measured at the same 
magnetic field value and plotted versus rig around the filling factor u = 2. The 
geometric capacitance [25] (the inverse of the first term in Eq. (2)), depicted 
by dashed line in Fig. 7(a), slightly increases with Ug since the thickness of the 
2D electron layer — and, therefore, the average distance between the 2D layer 
and the gate — decreases with the gate voltage. The second term in Eq. (2) 
is responsible for the dip centered at = 3.87 x 10^^ cm~^, corresponding to 
z/ = 2, and sharp maxima on both sides of it. Note that at these maxima, the 
capacitance exceeds Cq, which corresponds to the negative thermodynamic 
compressibility discussed above. Farther from integer filling factors, the intra- 
level interaction corrections become weak, being proportional to {z/}~^/^, and 
the measured capacitance approaches Cq (as long as the broadening of Landau 
levels is negligible, i.e., dus/dfi ^ dng/ dii\^^Q). 

Similar maxima on both sides of z/ = 2 are seen in the magnetization data 
shown in Fig. 7 (b). At the maxima, the derivative d^/dB exceeds maximum 
values possible in a non-interacting 2D electron gas, which are determined by 
Eq. (3) and are depicted in the figure by dashed lines. The possibility that 
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Fig. 8. Illustration of how the effective (7-factor (a, b) and the cyclotron mass (c, 
d) have been measured. The imaginary current component is plotted as a function 
of the deviation of the filling factor from = 2. In perpendicular magnetic fields, 
the difference between d^/dB for spin-down (|) and spin-up (j) electrons yields 
g* in units of the Bohr magneton. In tilted magnetic fields, the difference between 
dji/dB for electrons with = 1 and = is equal to 2iib (me /m*) cos cf). The 
dashed lines show noise-averaged values. 5B = 0.022 T (a, b) and 0.0055 T (c, d). 



dfi/dB might exceed its maximum non- interacting values due to intra- level 
Coulomb interactions between quasiparticles was predicted by MacDonald et 
al. [21]; in fact, this is how negative compressibility [22] manifests itself in 
magnetization measurements. Sharp spike just above z/ = 4 and maxima on 
both sides oi v = 2 m the dependence shown in Fig. 6 are of the same nature. 

It is straightforward to obtain the effective gf-factor from the data for dfi/dB. 
In accordance with Eq. (3), it is equal (in units of the Bohr magneton) to the 
difference between d^/dB for spin-down (|) and spin-up (t) electrons belong- 
ing to the same Landau level: fisg* = (dfi/dB) ^ — {dfi/dB)p It is important 
that this method of determining the (7-factor does not require the use of any 
fitting procedures or parameters. Figure 8 (a, b) shows measured d^/dB as a 
function of the deviation of the filling factor from 2 at two values of magnetic 
field. Near z/ = 2, there are sharp intra-level interaction-induced structures dis- 
cussed above; these regions have been excluded from the analysis. However, 
farther from u = 2, the dependences for u < 2 and u > 2 become parallel to 
each other. This ensures that the so-determined g* is not affected by the valley 
splitting [26,27] and intra-level interaction effects [21,22] discussed above. The 
latter contribute equally to both spin-up and spin-down dependences and can- 



12 



8 



7 



o 



6 



T3 
CO 



5 





m7m 



g72 



c 







0.5 



1.5 2.5 



3.5 



n (10" cm"^) 



Fig. 9. The effective (^-factor (circles) and the cyclotron mass (squares) as a func- 
tion of the electron density. The solid and long-dashed lines represent, respectively, 
the g(-factor and effective mass, previously obtained from transport measurements 
[28], and the dotted line is the Pauli spin susceptibility obtained by magnetization 
measurements in parallel magnetic fields. The critical density ric for the B = 
metal-insulator transition is indicated. 

eel eaeh other out. Disorder also contributes equally to dfi/dB on both sides 
of z/ = 2: we have found that at magnetic fields down to approximately 3 T, 
there are wide regions of filling factors where capacitance {i.e., the density of 
states) is symmetric around u = 2 (see, e.g., Fig. 7 (a)); furthermore, close- 
ness of the capacitance to Cq attests that the disorder- induced corrections are 
small. At lower magnetic fields, however, the electron-hole symmetry around 
z/ = 2 breaks down, which sets the lower boundary for the range of magnetic 
fields (and, consequently, electron densities). Note that temperature smears 
out the dependences in a way similar to disorder: at higher temperatures, the 
capacitance at half-integer filling factors decreases, which leads to a worsening 
of the method accuracy. 

In Fig. 9 we plot the measured (^-factor along with the one previously obtained 
from transport measurements (solid line). One can see that there is no sys- 
tematic dependence of the (^-factor on n^: it remains approximately constant 
and close to its value in bulk silicon even at the lowest electron densities, 
which is in good agreement with the transport [28] and magnetocapacitance 
[26] results. 

The same method can be used for determination of the cyclotron mass in tilted 
magnetic fields strong enough to completely polarize the electron spins [29]. 
If (and only if) the spin splitting exceeds the cyclotron splitting, the gap at 
z/ = 2 lies between Landau levels 0| and Ij, and the difference {dfj,/dB)]^=i — 
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{dii/dB)N=o is equal to 2fiB {iTie/m*) cos (p, where is the tilt angle. Once 
the electron spins are fully polarized at filling factors above z/ = 2, the tilt 
angle is automatically large enough for the level crossing to have occurred. The 
region of explorable electron densities is restricted from above by the condition 
that the electrons must be fully spin-polarized, while with our current set-up, 
the maximum magnetic field at which we can apply the modulation is only 
8 tesla capable of polarizing the electron spins up to ^ 2 x 10^^ cm^^ (see 
previous subsection and Ref. [30]). Figure 8 (c, d) shows dfj,/dB as a function 
of — 2| under the condition Ug < n* at two tilt angles [31]. The extracted 
cyclotron mass at electron densities 1.55 and 1.35 x 10^^ cm~^ is significantly 
enhanced. At densities below 1.35 x 10^^ cm~^, the symmetry of capacitance 
on both sides of the v = 2 gap breaks down, making the determination of 
m* impossible. As a result, we were only able to obtain two data points. 
Nevertheless, good agreement with the effective mass previously obtained by 
transport measurements (Fig. 9) demonstrates the applicability of the new 
method and adds credibility to both transport and magnetization results. 

We stress once again that the advantage of the new method we use here is 
that it allows determination of the spectrum of the 2D electron system under 
the condition that the Fermi level lies outside the spectral gaps, and the inter- 
level interactions are avoided. Being symmetric about z/ = 2, the intra-level 
interactions are canceled out in the data analysis and do not influence the 
extracted (^-factor and cyclotron mass. Therefore, the obtained values g* and 
m* are likely to be identical with those for a continuous spectrum, and the 
comparison with previously found values of the g^-factor and the effective mass 
is valid. 



4 Conclusions 

The Pauli spin susceptibility has been determined by measurements of the 
thermodynamic magnetization and density of states in a low-disordered, strongly 
correlated 2D electron system in silicon in parallel magnetic fields. It is found 
to behave critically near the zero-field MIT, which is characteristic of the exis- 
tence of a phase transition. Magnetization measurements in perpendicular and 
tilted magnetic fields allow determination of the spectrum characteristics of 2D 
electron systems and show that enhancement of the (^-factor is weak and prac- 
tically independent of the electron density, while the cyclotron mass becomes 
strongly enhanced as the density is decreased. The obtained data agree well 
with the gi-factor and effective mass obtained by transport measurements, as 
well as with the Pauli spin susceptibility obtained by magnetization measure- 
ments in parallel magnetic fields, even though the lowest electron densities 
reached in the second experiment are somewhat higher. Thus, we arrive at 
the conclusion that, unlike in the Stoner scenario, it is indeed the effective 
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mass that is responsible for the dramatically enhanced spin susceptibility at 
low electron densities. This result is consistent with conclusions of the re- 
cent theory by Punnoose and Finkelstein [32] who made a renormalization 
group analysis for multi- valley 2D systems and found that the effective mass 
dramatically increases at the metal-insulator transition while the ^^-factor re- 
mains nearly intact. We note, however, that the theoretical spin susceptibility 
diverges at disorder-dependent density [32], whereas the experimental x 
grows critically near the disorder-independent density [6]. 
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